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Abstract 

Mesophyll conductance (g m ) is now recognized as an important limiting process for photosynthesis, as it results in 
a significant decrease of C0 2 diffusion from substomatal cavities where water evaporation occurs, to chloroplast 
stroma. Over the past decade, an increasing number of studies proposed that g m can vary in the short term (e.g. 
minutes), but these variations are still controversial, especially those potentially induced by changing C0 2 and 
irradiance. In this study, g m data estimated with online 13 C discrimination recorded with a tunable diode laser 
absorption spectrometer (TDL-AS) during leaf gas exchange measurements, and based on the single point method, 
are presented. The data were obtained with three Eucalyptus species. A 50% decrease in g m was observed when the 
C0 2 mole fraction was increased from 300 ixmol moP 1 to 900 jjimol moP 1 , and a 60% increase when irradiance was 
increased from 200 fx.mo\ moP 1 to 1100 jjimol moP 1 photosynthetic photon flux density (PPFD). The relative 
contribution of respiration and photorespiration to overall 13 C discrimination was also estimated. Not taking this 
contribution into account may lead to a 50% underestimation of g m but had little effect on the C0 2 - and irradiance- 
induced changes. In conclusion, (i) the observed responses of g m to C0 2 and irradiance were not artefactual; (ii) the 
respiratory term is important to assess absolute values of gr m but has no impact on the responses to C0 2 and PPFD; 
and (iii) increasing irradiance and reducing the C0 2 mole fraction results in rapid increases in g m in Eucalyptus 
seedlings. 

Key words: Eucalyptus globulus, Eucalyptus saligna, Eucalyptus sieberii, photosynthesis, stomatal conductance, tunable diode 
laser absorption spectrometry. 



Introduction 

During photosynthesis, C0 2 diffuses from the atmosphere 
(at a mole fraction C a ) to the sites of carboxylation (C c ) 
inside the chloroplasts (Gaastra, 1959; Farquhar et al., 
1980). C0 2 crosses the leaf boundary layer and traverses 
stomatal pores into the substomatal cavity. C0 2 then 
diffuses through the gas phase between mesophyll cells 
before reaching cell walls, where it is solubilized. In the 



liquid phase, C0 2 crosses the plasma membrane, the 
cytosol, and finally the chloroplast membranes to reach 
the sites of carboxylation in the chloroplast stroma. Most of 
the carboxylation is by RubisCO (ribulose-l,5-bisphosphate 
carboxylase/oxygenase), while a small fraction (estimated at 
5%) is by phosphoewo/pyruvate carboxylase (PEPc) in the 
cytosol. The mesophyll conductance to C0 2 (g m ) represents 
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the conductance from the substomatal cavities at a mole 
fraction Q, to the sites of carboxylation, and includes gas 
and liquid phase transfer. 

In early leaf photosynthesis models, g m was considered to 
be infinite. This was an explicit assumption in the original 
formulation of Farquhar's model of C 3 photosynthesis. 
However, many subsequent studies showed that there was 
a difference between Q and C c , suggesting that g m might be 
finite (Evans et al, 1986; Lloyd et al, 1992) and that it 
affects estimation and interpretation of photoynthetic 
parameters such as the maximal carboxylation activity of 
RubisCO (F cmax ) or the maximal light-driven electron flux 
(4J (Epron et al, 1995; Niinemets et al, 2009a). g m can 
limit the rate of photosynthesis by 25-30% (Epron et al, 
1995) and this limitation may be as large as the one due to 
stomatal conductance (Warren and Adams, 2006; Flexas 
et al, 2007c). The range of g m variation among species is 
similar to that of stomatal conductance, from 0.005 mol 
m~ 2 s -1 up to 0.5 mol m~ 2 s -1 (see the reviews by Evans 
and Von Caemmerer, 1996; Flexas et al, 2008), with high 
g m occurring in herbaceous annuals and lower values in 
evergreen gymnosperms. g m is apparently influenced by leaf 
traits such as thickness, tissue density (and therefore leaf 
dry mass per unit area, LMA; for meta analyses, see Flexas 
et al, 2008; Niinemets et al, 2009/)), cell wall thickness, or 
the proportion of gaseous versus liquid mesophyll phases 
(Piel et al, 2002; Terashima et al, 2006; Evans et al, 2009). 
Such morphological leaf properties would confer stable g m 
in the short term (e.g. a few days). 

Since the early 1990s, g m was shown to be affected by 
environmental factors. At first, changes were thought to be 
rather long-term changes such as, for instance, the decline 
induced by leaf ageing (Scartazza et al, 1998), by a gradual 
drought stress (Comic et al, 1989; Roupsard et al, 1996; 
Warren, 20086), or by salinity (Bongi and Loreto, 1989). 
More recently, short-term changes were shown to occur in 
response to temperature (Bernacchi et al, 2002; Pons and 
Welschen, 2003; Warren and Dreyer, 2006; Warren, 2008a). 
Similarly rapid responses of g m (at the minutes scale) were 
reported to occur under varying C0 2 mole fractions. 
A review by Flexas et al (2008) reported that g m decreases 
with increasing C0 2 (data from Harley et al, 1992; During, 
2003; Flexas et al, 20076). This pattern was also observed 
by Hassiotou et al (2009) among seven Banksia species, by 
Vrabl et al. (2009) in Helianthus annus, and by Bunce (2010) 
in Glycine max and Phaseolus vulgaris. Tazoe et al. (2011) 
recently observed the same decreasing pattern among three 
species with contrasting photosynthetic capacities. Only 
a few studies have tested the rapid response of g m to 
irradiance. Flexas et al. (2008) reported an increase of g m 
with increasing irradiance (data from Gorton et al, 2003; 
Flexas et al, 20076). 

Nevertheless, there is still no consensus about the reality 
of such rapid responses of g m to C0 2 and irradiance. 
Several studies reported g m to be stable in response to 
changes in C0 2 (Von Caemmerer and Evans, 1991; Loreto 
et al, 1992; Tazoe et al, 2009) and irradiance (Tazoe et al, 
2009; Yamori et al, 2010). Some of the discrepancy may be 



due to measurement accuracy and/or artefacts. This can be 
true with the two main methods used to estimate g m : the 
fluorescence/gas exchange technique (Loreto et al, 1992) 
and the isotopic discrimination method (Pons et al, 2009). 
One of the complications with the isotopic method, which is 
used here, is that the contribution of 13 C discrimination 
during respiration and photorespiration needs to be taken 
into account (see the model description below). Some 
studies have ignored discrimination during respiration and 
photorespiration (Flexas et al, 20076, c; Vrabl et al, 2009), 
or approximated respiration in the light by that in the dark 
(Tazoe et al, 2009, 2011). The fractionation factors 
associated with respiration and photorespiration can be 
taken into consideration using recent estimates (Lanigan 
et al, 2008). An alternative approach is to limit the impact 
of fractionation during photorespiration by using low 0 2 in 
the measurement atmosphere (Tazoe et al, 2009). 

To date, there is no consensus regarding whether g m 
responds rapidly to irradiance and the C0 2 mole fraction, 
despite the importance of this for interpreting responses of 
photosynthesis to environmental variables. The aim of this 
research was to assess whether short-term variations of g m 
occur in response to changes in irradiance and the C0 2 
mole fraction. Rapid response of g m was assessed by 
recording online 13 C0 2 discrimination during photosynthe- 
sis with a custom-built photosynthesis chamber adapted to 
a LI-COR 6400, and coupled to a TDL-AS (or tunable 
diode laser absorption spectrometer). To examine the 
ubiquity of responses, three species of Eucalyptus with 
differing rates of photosynthesis were used, and the re- 
spiratory term of online 13 C discrimination was taken into 
consideration to reinforce the findings. Measurements were 
carried out on species with contrasting photosynthetic 
capacity to check whether g m variability is species de- 
pendent and how g m is related to A and g s within and 
among species. 

Materials and methods 

Plant material 

Three species in the genus Eucalyptus were used: E. globulus 
Labill., E. saligna Sm., and E. sieberi L. A. S. Johnson. Plants were 
grown for 4 months in a naturally illuminated greenhouse at an 
average daily temperature of 25 °C in 8.0 1 pots filled with 
compost-based substrate. They were watered with automatic drip 
irrigation. At the time of the experiment they were ~60 cm high. 

Gas exchange measurements 

The response of the net C0 2 assimilation rate (A, umol m~ 2 s _1 ) to 
variations in substomatal CO2 mole fraction (Q) and photosyn- 
thetic photon flux density (PPFD, umol mT 2 s ) was measured in 
each of four replicate plants of the three species, for each 
treatment (e.g. 2 treatmentsx3 speciesx4 replicates=24 plants). 
All measurements were made on the youngest fully expanded 
leaves. Plants were transferred from the greenhouse to the 
laboratory at 25 °C. To assess photosynthesis and discrimination 
against 1 C0 2 simultaneously, a LI-6400 portable gas exchange 
system (LI-COR, Lincoln, NE, USA) equipped with a custom- 
built chamber of 18 cm 2 and coupled to a TDL-AS (TGA100A, 
Campbell Scientific, Logan, UT, USA) was used. 



Mesophyll conductance varies with C0 2 and irradiance | 5337 



Before measurements, leaves were exposed to 1000 umol rrT 2 s _1 
PPFD and 400 umol C0 2 moL 1 for 30 min to induce photosynthe- 
sis and stomatal opening. C0 2 responses of photosynthesis were 
measured under a constant PPFD of 1000 umol irT 2 s~ while the 
C0 2 mole fraction for the reference gas (C e ) was sequentially set at 
300, 500, 700, and 900 umol mol" 1 . The PPFD responses at 200, 
500, 800, and 1100 umol rrT 2 s were measured under a constant 
C0 2 mole fraction of 400 umol mol - . During each step of the C0 2 
or PPFD responses, at least 20 min was allowed for A and g s to 
reach a steady state, and three individual points separated by 180 s 
were recorded. Each set of three individual data points was 
averaged to characterize the response at one step in the response 
curve. The standard error of the three individual points was used to 
assess the analytical variability. 

Air flow into the chamber was set at 400 umol s _1 , and leaf 
temperature at 25 °C. Irradiance was provided by a LI-COR RGB 
light source (6400-18, LI-COR-0128), which covered the entire 
chamber surface. A subsample of air from the sample and 
reference gas lines of the LI-6400 was diverted to the TDL-AS. 
The TDL-AS measured sequentially gas from two calibration 
tanks, the LI-6400 reference gas, and finally the LI-6400 sample 
gas. Each intake was measured for 45 s, with the first 15 s ignored 
to minimize carryover and enable stabilization between intakes. 
The total time for a sequence of measurements was therefore 180 s. 

The TDL-AS gas was connected through a 'T' tubing to the 
reference tube of the LI-6400 between the console and the IRGA. 
In the same way, the TDL-AS intake of sample gas was connected 
in between the LI-6400 chamber exhaust. The TDL-AS was set to 
continuously withdraw 150 ml min -1 ( — 102 umol s _1 ) from each 
of the sample and reference fluxes of the LI-6400. This withdrawal 
of air was much smaller than the flow through the LI-6400 
chamber, which means that the TDL-AS could sample air from 
the LI-6400 while maintaining a positive pressure in the LI-6400 
chamber. The existence of a positive pressure inside the LI-6400 
chamber was checked through the curvature of the propafilm 
covering the top of the chamber. TDL-AS data were matched with 
LI-6400 data by taking into account the time lag of 37 s that was 
recorded between the chamber and the TDL-AS. The TDL-AS 
data were only used to record the isotopic composition of the 
reference and sample gas, while all photosynthesis parameters were 
estimated from the LI-6400 data. 



Isotopic measurements and system testing 

Discrimination by a leaf was assessed by measuring the isotopic 
composition in reference (S 13 C e , % 0 ) and sample gas (8 13 C 0 , % 0 ). 
The isotopic composition (8 13 C) was expressed as: 



Rs 



1 X1000 



where R s is the isotopic ratio (i?= 13 C0 2 / 12 C0 2 ) of the sample 
and ^vpdb is the isotopic ratio of Vienna Pee Dee Belemnite 
(VPDB, 0.0112372). 

The TDL-AS was calibrated with two tanks (Tl and T2) with 
C0 2 concentrations of 419± 10 umol moL 1 and 290 ±7 umol 
moP 1 [mean ± the confidence interval (CI)], respectively, given by 
the provider and checked with a recently factory-calibrated LI- 
8100 IRGA (LI-COR). The isotopic composition of the each tank 
was measured by sampling air into 12 ml exetainers (Labco 
Limited, Buckinghamshire, UK) and analysed via the gas-bench 
inlet of an IRMS (Delta S, Finningan, Bremen, Germany) at 
INRA Nancy (w=10 exetainers per tank). For Tl and T2, 5 13 C 
was -36.6±0.08% o and -36.9±0.2% o [mean ± standard deviation 
(Sd), »=10], respectively. Absolute values of 12 CO z and 13 C0 2 
were, respectively, 414.51 umol moL 1 and 4.487 umol moF 1 for 
Tl and 286.89 umol moL 1 and 3.104 umol moL 1 for T2, 
considering the C0 2 mole fraction indicated by the provider took 
into consideration both isotopologues. A linear interpolation was 



used between these two points for each isotopologue. For further 
calibration, a range of C0 2 mole fractions of 200, 300, 500, 700, 
1000, 1500, and 2000 umol moL 1 was generated using the C0 2 
injector of the LI-6400 fed with the same C0 2 cartridge during the 
whole test. C0 2 mole fractions above the calibration range led to 
a small deviation of +2.5% 0 . A second order polynom was fitted 
to describe the deviation of apparent 5 13 C from reference values 
measured at C e =300 umol moP 1 (because this is within 
the calibration range) along the extended C0 2 range 
(8 13 C=0.000007[CO 2 ] 2 -0.0035[CO 2 ]-5.2922; R 2 = 0.92, «=140). 
All TDL-AS values were corrected and a stable 8 13 C signal was 
obtained along the extended C0 2 mole fraction range. 

The noise of the system was assessed by observing the Sd of 
8 13 C values within each C0 2 step. An average Sd of 0.16 % 0 was 
observed for 8 13 C («=21). These values were used for the 
computation of the standard deviation of the observed discrimina- 
tion by the leaf (Sd Ao bs> see below). 

Finally the empty photosynthesis chamber was used to test the 
absence of a 8 13 C difference between the inlet and outlet that might 
be caused by leaks and/or C0 2 adsorption/desorption processes. It 
was observed that the 8 13 C e -8 13 C 0 difference was stable along the 
C0 2 mole fraction range used during the test, and it was confirmed 
that the empty chamber did not affect the isotopic composition of 
the air (i.e. leaks and C0 2 adsorption/desorption processes were 
negligible). The observed 8 13 C e -S C Q difference was much smaller 
than the 8 13 C difference recorded during actual measurements 
(Fig. 1), although under specific conditions (low photosynthesis) the 
two could overlap. A data filtering procedure was therefore 
implemented (see below). 
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Fig. 1. Boxplots of the difference in 8 C between the outlet and 
inlet of the 1 8 cm 2 chamber for (i) an empty chamber with inlet air 
C0 2 varying from 200 umol mor 1 to 1000 umol mor 1 (n=10); (ii) 
the whole set of measurements under varying C0 2 concentrations 
(300 umol moL 1 to 900 umol moL 1 ), n=40; and (iii) the whole set 
of measurements under varying PPFD (from 200 umol rrT 2 s~ 1 up 
to 1 100 umol rrT 2 s~ 1 ), n=40. The middle line represents the 
median, the upper and lower box limit the 75% and 25% quartiles, 
respectively, and whiskers represent the extreme value. 
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Model description 

The observed discrimination (A obs ) is usually calculated following 
Evans et al. (1986): 



*obs - 



t(6 u C 0 - 



<5 13 c; 



1000 + <5 13 C o 



a{s n c 0 



<5 13 C, 



where: 



(i) 



(2) 



% is the ratio of C0 2 entering the chamber over the C0 2 
drawdown induced by the leaf. 

A is the result of discrimination by diffusion processes during 
C0 2 movement from the atmosphere to the chloroplast, and 
biochemical fractionation during carboxylation processes. Each 
fractionation step is characterized by a fractionation factor (due to 
diffusion or biochemistry) weighted by the gradient of concentra- 
tion. In the complete form, A is predicted by (Evans et al, 1986): 



C'Ci / i 7 

^ = a h — h a — h \e s + a, — h ft— 



(3) 



where: 



• a h is the fractionation during C0 2 diffusion in the boundary 
layer (2.9%,, Evans et al, 1986); 

• a is the fractionation during C0 2 diffusion in air through 
stomata into the leaf (4.4% 0 , O'Leary, 1981); 

• e s is the fractionation occurring when C0 2 is dissolved in the 
cell solution (1.1% 0 at 25 °C, O'Leary, 1981); 

• flj is the fractionation during C0 2 diffusion in the liquid phase 
(0.7% o , O'Leary, 1981); 

• b is the discrimination during carboxylation, and is dependent 
on fractionation by both RubisCO (b 3 =30°/ oo ) and PEPc (b 4 = 
-5.7% 0 ), 

b is computed as: 



b = (1 - j8)fej - /» 4 



(4) 



where (3 (between 0.05 and 0.1) is the relative amount of carbon 
fixed by PEPc (Farquhar and Richards, 1984). In the present 
experiment a value of Z>=28% 0 was used; that is, a P value of 
0.055. 

• / denotes overall discrimination during photorespiration. The 
value of / was set at 1 1% 0 , according to the theoretical approach 
of Tcherkez (2006), which was confirmed later by Lanigan et al. 
(2008); 

• e denotes overall fractionation during day respiration relative 
to photosynthetic products (Rd), and can vary between -10% o 
and +10% o (Ghashghaie et al, 2003). e was set at \% before 
correction following Wingate et al. (2007); 

• k is the carboxylation efficiency computed as (Farquhar et al. , 
1982); 



decarboxylation terms are ignored and g m is considered to be 
infinite: 



Aj = a+ b — a 



C a 



(6) 



With this 'single point method' first developed by Lloyd et al. 
(1992) and recently described by Pons et al (2009), g m can be 
estimated from a single value of A obs : 



(b - e s - a i)A/C a 

eRj/k+fr* 



(At - A ob )j 



(7) 



Model parameters 

Rd and T* were estimated with the 'Laisk method' (Viil et al, 
1977) for the three species (i.e. from the intersection of three 
A-Ci curves recorded at PPFDs of 100, 50, and 25 umol photon 
mT 2 s _1 and C e of 125, 100, and 50 umol moL 1 ). The 'Laisk 
method' provides Q* or the 'apparent' C0 2 compensation point in 
the absence of day respiration (Von Caemmerer et al, 1994), and 
was used as a proxy of r*. Because these measurements are 
sensitive to errors due to C0 2 leak diffusion (low A and low C a 
compared with the atmosphere), the potential C0 2 leaks due to 
diffusion though chamber gaskets were estimated (Flexas et al, 
2007a; Rodeghiero et al, 2007). A diffusion coefficient of the 
gaskets was computed with the procedure provided in the user 
manual of LI-COR, and a value of 0.938 umol s~ was found 
(while it usually is 0.46 for smaller 6 cm 2 chambers). This 
correction was incorporated into all gas exchange computations 
used for Y* and Rd estimations. The computed values of T* did 
not differ between species, so the mean (r*=38.7±0.51 umol 
mol -1 , w=13) was used as a common value for the three species. 
For E. globulus, E. saligna, and E. sieberi, R A was 0.41 ±0.09, 
0.31±0.09, and 0.68±0.07 umol m~ 2 s~\ respectively. 

As the isotopic signature of the reference gas provided by 
cartridges differed from that used by the leaves for earlier 
photosynthesis, e was replaced by e'=e+8 13 C tan k-5 13 C at mosphere 
(Wingate et al, 2007). 5 13 C in the cartridge was measured with the 
TDL-AS at the chamber inlet. In the present case, the LI-6400 was 
fed with compressed C0 2 cartridges with 5 13 C tank varying between 
-1% 0 and -4% 0 except for two cartridges with 5 13 C tan k= -19% 0 . <?' 
therefore varied between +4% 0 and +6% 0 , except for two plants 
where it was -11% 0 . 

C c , the C0 2 concentration at the site of carboxylation, was 
calculated from Fick's Law as: 



Q- 



(8) 



Finally, the total leaf conductance to C0 2 was calculated 
following Ball (1988), assuming resistances in series as: 



1 



where the stomatal conductance to C0 2 is g sc =g s Jl.6. 



(9) 



k = (A+R^KC, - r) (5) 



• T* is the C0 2 compensation point in the absence of day 
respiration (Brooks and Farquhar, 1985). 

The estimation of g m is based on the difference between the 
observed discrimination by the leaf (A obs ) and the discrimination 
predicted from the simplified form of the model (A;) in which 



Propagation of uncertainty from measurement to A calculation and 
data filtering 

The uncertainty (standard deviation) of A obs due to the finite 
precision of 8 13 C measurements was estimated. This was achieved 
by propagating uncertainty (standard deviations) of 8 13 C e and 
5 3 C Q through the equations estimating A obs (see Appendix for 
details): 
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AObS ~{ W»C 0 -6»C.) (iq) 

1 + V Sdsnco 2 — Sdsnce 2 + Sdsnco* ] ^ 
l+S u C 0 -H(S u C 0 -S u Q) J 

The propagated uncertainty in A obs (i.e. SdAobs) was used as the 
basis for a filter to remove unreliable estimates of A obs . Computa- 
tion of gm is based on the difference between A; and A obs (i.e. 
Ai-A obs ), thus it was reasoned that g m estimates would be unreli- 
able if the difference A ; -A obs was smaller than Sd Ao bs- Conse- 
quently all values where A obs +Sd Aobs >A; were rejected. This filter 
was applied to the individual points in the data set, rejecting 33 
among 238 points. 

Statistical analyses 

All statistical analyses were performed with R (R Development 
Core Team 2010, http://www.R-project.org). Mixed-effect linear 
models were run to assess species and treatment effects on A, g s , 
g m , A;-A, C c , and the Q-C c drawdown, as shown in Table 1. For 
the C0 2 treatment, 'species' (as a factor) and 'Q' (as a covariate) 
were incorporated into the model as fixed effects, and 'individual 
within species' as a random effect. For variations of C c and the 
Q-C c drawdown, C a was used as covariate. For PPFD treatment, 
species and PPFD were set as factors. Normality and heterosce- 
dasticity were graphically checked with QQ-plots. In the case of 
heteroscedastic data, the mean was weighted as a function of the 
variance. In the case of non-normal distribution, variables were 
log-transformed. The species x treatment interaction was tested for 
each procedure, and was removed from the model when not 
significant. In the absence of interaction, comparison of the 
intercepts was performed to assess differences between species. In 
the case of interaction, slope comparisons were performed to test if 
species responses differed from each other. Significance was 
accepted at P <0.05. Mean least squares regression was used to 
assess the correlation between variables (R 2 and P-value). 



Results 

Variation of g m under changing C e 

The C0 2 mole fraction was changed in the air entering 
the chamber in three steps from 300 umol mol -1 to 900 
umol mol -1 , inducing a range of Q from 185 umol mol -1 to 
745 umol mol -1 . Net C0 2 assimilation rate (^4) was 
positively related to Q and varied between 3 umol mol -2 
s -1 and 18 umol mol -2 s -1 , while stomatal conductance to 
water vapour (g s ) was negatively related to Q and varied 



between 0.02 mol m" 2 s" 1 and 0.8 mol m" 2 s" 1 (Fig. 2). 
Eucalyptus sieberi had significantly higher A and g s than 
E. globulus and E. saligna (/-test P <0.05, Fig. 2). 

The difference between 8 13 C e and § 13 C 0 decreased with 
increasing Q (data not shown). Among all species, A obs 
varied between 12% c and 22% c and was positively correlated 
to Q/C a (7? 2 =0.79, P <0.001, data not shown). The 
difference between A calculated with infinite g m and no 
respiratory term (simple model) and observed A (A ; -A obs ) 
varied between 3%„ and 7% 0 but showed no clear trend with 
Q (Fig. 2, Table 1). 

g m computed by taking into account the respiratory 
component of discrimination varied from 0.025 mol 
m~ 2 s _1 to 0.55 mol m~ 2 s -1 (Fig. 2). g m was larger in 
E. sieberi than in the other two species (r-test, P <0.05). g m 
was affected by Q (Table 1), and decreased when Q 
increased. Post-hoc tests revealed that E. globulus and 
E. saligna displayed g m -Q slopes significantly different 
from zero. In E. sieberi, three individuals out of four 
showed a clear negative pattern when Q increased, but not 
the fourth. The relationship between g m and g s was 
significant among all species (R 2 ~0.54, P <0.001), and 
within E. globulus and E. saligna when treated separately 
(data not shown). 

C c and the Q-C c drawdown were, as expected, severely 
affected by C a (Table 1). The Q-C c drawdown was —50 
umol mol -1 at C a =200 umol mol -1 and increased up to 200 
umol mol -1 at C a = 900 umol mol -1 (Fig. 3). When tested 
individually, all slopes of the responses of C c and Q-C c to 
Q were different from zero. 

Variation of g m with changing irradiance 

A and g s increased significantly with irradiance (Table 1). 
A and g s were larger in E. sieberi than in the two other 
species (see Fig. 4, f-test P <0.05). In each species, g s 
increased significantly from 200 umol m~ 2 s -1 to 500 umol 
m~ 2 s -1 PPFD, and then stabilized (Fig. 4). A and g s were 
positively correlated, for all species taken together (Fig. 5C, 
R 2 =0.16, P <0.001). Across the range of irradiance, A obs 
varied between 1 1% C and 20% o and was positively correlated 
with Q/C a CR 2 =0.74, P <0.001, data not shown). A ; -A obs 
varied between 3% 0 and 9%„ but was not affected by 
irradiance or by species (Table 1, Fig. 4). 



Table 1. Mixed effects model for A, g s , g m , and A r A 

Species, C a and C,, or PPFD effects were incorporated into the model as fixed effects, and individual plant as a random effect. In the case of 
heteroscedastic data the mean was weighted as a function of the variance. For Q and C a , the degree of freedom (df) was 1 , for PPFD df=3, 
and for species df=2. Significant values (P <0.05) are shown in bold. 
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NS 
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NS 


0.008 


0.038 
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NS 
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Fig. 2. Relationships between the C0 2 mole fraction in the substomatal cavities (Q) and (A) net C0 2 assimilation rate (A). (B) Stomatal 
conductance to water vapour (g s ). (C) Difference between predicted and measured isotopic discrimination (A-A). (D) Mesophyll 
conductance (g m ). Eucalyptus globulus is in black, £ saligna in grey, and E. sieberi in white. The SE is provided by the average of three 
measurements taken at 180 s intervals. Measurements were made at four levels of C e on four plants per species, and were filtered 
against noisy values of 5 13 C e -5 13 C 0 . 
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Fig. 3. Left: C0 2 mole fraction in the substomatal cavities (Q, disks) and at carboxylation sites (C c , triangles) as a function of C a . The 
dashed grey line is the 1:1 relationship. Right: C r C c drawdown as a function of C a . Each point represents the mean of 1-3 analytical 
measurements ±SE, with £ globulus in black, £ saligna in grey, and £ sieberi in white. 



g m varied between 0.04 mol mol 2 s 1 and 0.6 mol mol 2 
s -1 , and was positively related to PPFD. As for g s , the 



response consisted of a significant increase between 200 



umol 



m 



s and 500 umol 



m 



PPFD with 
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Fig. 4. (A) Net C0 2 assimilation rate (4). (B) Stomatal conductance to water vapour (g g ). (C) Difference between predicted and 
measured isotopic discrimination (A r A). (D) Mesophyll conductance (g m ). (E) C0 2 mole fraction at carboxylation sites (C c ). (D) The CrC G 
drawdown at four different levels of PPFD. Means ±SE (n=2-4 replicate plants), with £ globulus in black, E. saligna in grey, and 
£ sieberi in white. 
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Fig. 5. Relationships between net C0 2 assimilation rate (/\) and (A) total leaf conductance to C0 2 (g t , black squares), (B) mesophyll 
conductance to C0 2 (g m , grey squares), and (C) stomatal conductance to C0 2 (g sc , white squares) under varying PPFD. 



a stabilization above this threshold. g m was positively 
correlated with g s (R 2 =036, P <0.00T, data not shown) 
and A (R 2 =0A9, P <0.001, Fig. 5B), among all species. 
Significant A-g m and g m -g s relationships within each species 
were also detected (except for g s -g m in E. saligna). 



Total leaf conductance to C0 2 (g t ) was strongly corre- 
lated to A CR 2 =0.83, P <0.00T) as shown in Fig. 5A. Over 
the full set of irradiance values, C c varied between 200 (imol 
mol -1 and 260 umol mol -1 , and the Q-C c drawdown 
between 40 umol mol -1 and 60 umol mol -1 . None of these 
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parameters displayed any variation with irradiance (Table 
1), or with species. 



Discussion 

This study provides support to recent evidence that g m 
varies rapidly (within minutes) in response to environmental 
conditions. The rapid responses were observed under two 
different sources of variation: C0 2 mole fraction and 
irradiance. In seedlings of three Eucalyptus species a modest 
but significant decrease of g m with increasing C0 2 mole 
fraction, and a significant increase with irradiance, was 
found. The effect was visible in the three species irrespective 
of the photosynthetic capacity. 

Importance of the respiratory and photorespiratory 
terms in the estimation of g m 

The isotopic method estimates g m from the 13 C0 2 discrim- 
ination during photosynthesis by comparing observed 
values with those derived from a model-based prediction of 
discrimination under infinite mesophyll conductance. This 
approach requires a high precision in discrimination 
records, which is now achieved by combining precise leaf 
gas exchange measurements with online TDL-AS records of 
changes in 13 C0 2 / 12 C0 2 in the atmosphere around the leaf 
(for a discussion of the technique, see Pons et al, 2009; 
Tazoe et al, 2011). One of the important problems with this 
method is the fact that several discrimination steps during 
photosynthesis, respiration, and photorespiration need be 
taken into account. In particular, a number of earlier 
studies omitted the respiratory and photorespiratory terms 
(Flexas et al, 2007b; Vrabl et al, 2009). Moreover, the 
response of g m to C0 2 was affected by the 0 2 mole fraction 
in the air; that is, by the occurrence of photorespiration 
during measurements: it was visible only under low 0 2 
(Tazoe et al, 201 1). 

In the present study, these two terms were incorporated 
into the g m estimates displayed in the results. Absolute 
values of g m were up to 50% larger when these terms were 
incorporated, and this enhancement was independent of the 
treatments applied (Fig. 6). Whether this potentially large 
change had an impact on the observed the effects of 
changing the C0 2 mole fraction in the air or irradiance was 
tested: the response of g m to C0 2 mole fraction and 
irradiance remained significant even when the respiratory 
term was omitted (see Table 2). the effects of substituting e 
(fractionation during day respiration) and / (fractionation 
during photorespiration) with extreme values (f—0% o , e= 
-10% o , e= +10% o ) were similarly tested. Despite the fact 
that these changes resulted in significant differences of 
computed g m , they did not result in any loss of significance 
of the observed effects of C0 2 or irradiance (data not 
shown). Vrabl et al. (2009) compared g m estimated with the 
fluorescence method (which takes the respiratory terms 
implicitly into account) and with the isotopic method 
(without taking them into account) and found the same 



range of g m values and the same negative response to Q. 
The present computations suggest that the respiratory terms 
can be important in the estimation of the absolute values of 
g m , but have only little influence on the observed C0 2 or 
irradiance responses of g m . The respiratory terms of isotopic 
discrimination are unlikely to be responsible for the 
discrepancies among studies. Addressing the question of 
changing e and / during the C0 2 and irradiance treatments 
was omitted in this discussion: up to now there is no reason 
to assume that these values are not stable across the whole 
range of environmental conditions. 



Table 2. Impact of omitting the respiratory term in the discrimina- 
tion model on the assessment of the impact of the C0 2 mole 
fraction (Q) or of irradiance on g m 

Mixed effect model for mesophyll conductance computed including 
(as shown in Table 1) or omitting the respiration and photorespiration 
terms. 
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9m 


g m (omitting the 
respiration terms) 


9m 


g m (omitting the 
respiration terms) 


Variable 
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18.67 


13.14 


6.41 


12.09 




P 


<0.001 


0.001 


0.002 


<0.001 


Species 


F 


10.11 


10.28 


NS 


4.89 
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0.006 


0.006 


NS 


0.03 


Interaction 
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3.44 


3.64 


NS 


NS 




P 


0.047 


0.04 


NS 


NS 




Mesophyll conductance without respiration term 
(mol nrV 1 ) 

Fig. 6. Mesophyll conductance g m as computed by taking into 
account the contribution of respiration and photorespiration to 
13 C0 2 discrimination versus without taking them into account (i.e. 
fractionation factors e and f set to 0 in Equation 7). Results 
obtained under changing C0 2 (filled squares) or irradiance (open 
squares). The dashed line represents the 1 :1 relationship. Each 
point represents the mean ±SE of a given measurement (n=3). 
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Response of g m to C0 2 mole fraction 

There was at maximum a 50% decrease in g m with 
increasing C0 2 mole fraction (from 300 umol mol -1 to 900 
umol mol -1 ). This response was general, as all three species 
displayed a similar response to C0 2 , with the exception of 
a single individual of E. sieberi. The results contrast with 
earlier studies reporting no response of g m to C0 2 . This was 
the case in Quercus ilex and Citrus aurantium (chlorophyll 
fluorescence method; Loreto et al, 1992), Raphanus sativus 
(isotopic discrimination; Von Caemmerer and Evans, 1991), 
and Triticum aestivum (isotopic discrimination; Tazoe et al. , 
2009). The present results confirm several studies that 
reported a negative relationship between g m and C0 2 
(Flexas et al, 2007ft; Hassiotou et al, 2009; Vrabl et al, 
2009; Bunce, 2010). The same magnitude of decrease of g m 
with Q as in Flexas et al. (2007ft) and Vrabl et al. (2009) 
was found. Tazoe et al. (2011) found an —30% decrease of 
g m when Q increased under 1% 0 2 but no effect under 21% 
0 2 in two Arabidopsis thaliana genotypes or in Nicotinia 
tabacum, while in the present study a significant effect was 
found even under 21% 0 2 . The absence of a clear pattern of 
g m responses among studies could be interpreted as 
a species-dependent response of g m to Q, but no common 
trait seems to be shared by the 'non-responsive' versus the 
'responsive' species. 

There are a suite of methodological reasons and 
artefacts (e.g. choice or calculation of ft, e, f, R d , or F*) 
that might explain the discrepancy among studies. The 
case for e and / was discussed above. Measured values of 
Rd and T* (i.e. of its proxy Q*) were used for each species 
rather than arbitrary values taken from the literature. 
Despite some uncertainties regarding the choice of ft or 
/ (Lanigan et al, 2008; Pons et al, 2009), it is concluded 
that the response of g m that was recorded here is unlikely 
to be an artefact. 

Rapid responses of g m to C0 2 could be mediated by 
aquaporins that might impact the permeability of plasma 
and chloroplast membranes to C0 2 (Terashima and Ono, 
2002), enhancing C0 2 diffusion in the liquid phase. Flexas 
et al. (2006) stated that the expression of NtAQPl 
aquaporins can change g m values by 20-50% in N. tabacum. 
These variations of g m are of the same magnitude as those 
in the present study, but only direct measurements of 
aquaporin expression/activity could confirm the role of 
these proteins in the diffusion of C0 2 through membranes. 
Establishing a parallel between short-term responses of g m 
and of the expression and activities of aquaporins is still an 
open area for research. Tholen et al. (2008) also found that 
chloroplast movements can induce a variation of g m by 50% 
in A. thaliana. However, it is not known yet whether C0 2 
variations can directly mediate a displacement of chloro- 
plasts. 

A positive relationship was observed between g m and 
g s , as was also observed by Flexas et al. (2007c). In- 
terestingly, such a relationship was also found by Flexas 
et al. (2006) when they compared plants overexpressing 
NtAQPl aquaporin and controls. They suggested that 



the variation of g m primarily induced by manipulating 
NtAQPl expression probably also led to an adjustment of 
g s and subsequently of A. This potentially indicates 
a physiological link between these two parameters (Flexas 
et al, 2007c; Vrabl et al, 2009). Nevertheless, the precise 
signalling cascade that could cause a coordinated response 
of stomatal and of mesophyll conductance remains to be 
elucidated. 

Response of g m to irradiance 

Several studies observed an increase in g m with increasing 
irradiance. Flexas et al. (2007ft) reported that g m increased 
in tobacco by —40% when irradiance increased from 250 
umol m~ 2 s _1 to 1000 umol m -2 s~ . Data from Gorton 
et al. (2003), reanalysed by Flexas et al. (2008), also showed 
a positive effect of irradiance on g m . Hassiotou et al. (2009) 
detected an —22% increase in six Banksia species as 
irradiance was switched from 500 umol m -2 s _1 to 1500 
umol m -2 s~ . The present study corroborates this positive 
effect of irradiance, with an increase in g m by 60%, up to 
a plateau in irradiance reached between 500 umol m~ 2 s -1 
and 1100 umol m -2 s~ , depending on the species. A slightly 
larger sensitivity of g m to irradiance than earlier studies was 
found. Sensitivity to irradiance could be (i) species de- 
pendent or (ii) due to different parameterization of the 
model enabling g m estimation. For instance, in the present 
case, removing the respiratory and photorespiratory terms 
in the estimation of g m led to a slightly smaller response of 
g m to irradiance (Table 2). A systematic analysis of the 
reported responses, with a standardized parameterization, 
would be very helpful. 

As during the response to C0 2 , g m was positively 
correlated to g s . This relationship seems to be independent 
of the method used to vary net C0 2 assimilation rates. 
A review by Flexas et al. (2008) insisted that g s and g m 
responded in parallel to irradiance, C0 2 , temperature, and 
drought stress (Warren, 2008ft). On the other hand, Warren 
(2008ft) reported that g m was unaffected by increases in 
vapour pressure deficit (VPD) while g s decreased strongly, 
and Vrabl et al. (2009) reported that g m was unaffected by 
feeding leaves with abscisic acid (ABA), whereas there was 
a clear decrease in g s . In the two latter cases the net C0 2 
assimilation rate remained unaffected by VPD and ABA 
despite the severe reduction of g s . The g s -g m relationship 
therefore may reflect a tight coordination between A and 
g m . Thus it seems that g m contributes to adjust the C0 2 
supply to the sites of carboxylation in response to photo- 
synthetic limitations such as light availability, hydraulic 
constraints, or biochemical limitations (Warren et al., 2007). 
In the present study, the coordinated variations of g s and g m 
apparently led to a very stable C c (and Cj-C c drawdown) 
across irradiance levels, despite large variation in A. Such 
a homeostasis of C c was already observed across a range of 
leaf morphologies (Hassiotou et al, 2009), during leaf 
ageing (Ethier et al, 2006; Warren, 2006), and such an 
adjustment seems also to occur during short-term fluctua- 
tions of irradiance. 
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Conclusion 

This study with three Eucalyptus species confirmed that g m 
estimated with the online 13 C discrimination method 
declines in response to short-term increases of the C0 2 
mole fraction and increases with irradiance. The response 
to irradiance is saturated above 500 umol m -2 s -1 PPFD. 
The respiratory term in the 13 C discrimination equation 
was found to be important to estimate absolute values of 
g m but had little impact on the C0 2 and irradiance 
responses. During the responses to C0 2 and PPFD, g s and 
g m were tightly correlated and varied in parallel indepen- 
dently of the source of variation. Moreover, it was 
observed that coordinated adjustments of C0 2 demand 
(A) and supply (g s and g m ) led to a stability of C c across 
irradiance variations. C c homeostasis could be an advan- 
tage for the leaf to prevent large variation of the 
oxygenation/carboxylation ratio of the RubisCO, when 
the C0 2 demand increases. 
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Appendix I 

Here we describe the procedure to calculate the standard 
deviation of A obs (SdAobs) based on the standard deviation of the 
isotopic composition of reference and measured air (Sdsi 3 c e and 
Sdsnco: respectively). In the equation of Evans et al. (1986) to 
calculate A obs , 8 13 C e and 8 13 C 0 were substituted by their respective 
standard errors Sd 513Ce and Sd 513Co to obtain: 



Sd, 



Aobs 



5UCo 



Sd. 



S\3Ce 



1 + Sdsnco - ^(Sdsnco - Sdsnce) 



(Sd 8 



o~~ Sdgisce and Sd 513Co were replaced by 



y'Sdsnco + Sdsnce and y/Sdma,, respectively, as is described in 
Harris (1991) to obtain: 
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1 3 Co 
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Then, to calculate the cumulative effect of the upper term and 
the lower term, their relative Sds were summed to obtain the 
relative standard error of A obs : 
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SdAobs is then equal to: 
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